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^^ Abstract 

o , 

^^ ■ It is well understood that the studies of correlations between produced par- 

^~>! tides, the effects of coherence and chaoticity, an estimation of particle emit- 

ting source size play an important role in high energy physics [1]. We mean 
the investigation of the space-time extension or even squeezing of particle 
sources via the multiparticle quantum-statistics correlation. We obtain the 
two-photon correlation function that can provide the space-time information 

r^ • about the Higgs-boson source in thermal environment and we argue that such 

an investigation could probe the Higgs-boson mass. 
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1. The Large Hadron Collider (LHC) at CERN is at the stage to provide particle 
^ I physicists with treasure of data. These data allowed a precise measurement of many 

ly-s ' important parameters of modern particle physics in order to test their consistency 

O ' and to discover the Higgs-boson. At the same time, Uttle is known about Higgs-gauge 

bosons interplay in particular in theoretical aspects. The effect of Bose-Einstein 
correlations (BEC) is clear and undeniable part of this theory, complicating the 
O I quantum statistical description of multi-lepton final states. 

Historically, the BEC measurements have concentrated on pion pairs correla- 
tions, however have also been applied to more heavy hadrons, quarks, protons, and 
^ ■ even gauge bosons, photons, Z-bosons etc. 

It is known that the wave function describing a system of two identical bosons 
should be symmetric under permutation of these bosons. As a consequence, the 
four-momentum differences between two bosons will be smaller than in the world 
where Bose-Einstein statistics would not apply. 

Actually, the direct photon BEC can provide information about the space-time 
distribution of the hot matter prior to freeze-out. However, the BEC with direct 
(primary) photons is faced with difficulties compared to hadron BEC primarily due 
to the small yield of photons emitted directly from the hot region just after particle 
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collisions. The main background of photons is produced by decays of final hadrons 
or gauge bosons, or even (pseudo) scalar particles. 

In Higgs searching program at the LHC, we suppose that two photons are pro- 
duced mainly through Higgs-boson decay. From the theoretical point of view the 
decay of Higgs into two photons emerges through the one-loop pattern: the Higgs- 
boson transferred in pair of a quark and antiquark or vector bosons, or even in scalar 
particles. 

In order to reject most of the hadron background, the reconstructed events in 
a spectrometer were required to have a deposit energy (level) of greater than the 
energy well above the minimal energy at which photons emerge from hadrons. Ex- 
perimentally, the vr" background constitutes a major difficulty which, however, to 
some extent, can be taken care of by measuring the photon pair invariant mass. 

The two-photon correlations in Higgs physics provide a powerful tool to explore 
the Higgs-boson mass estimated via the correlation radius which defines the geo- 
metrical size of the two-photon source. Hence, the two-photon correlation function 
is strongly dependent on the space-time properties of the Higgs decays. 

In the quark-loop scheme, the non-relativistic bremsstrahlung formula for the 
current in quark- antiquark interaction for one photon emission (with four-momentum 
k = (fcO, k)) is 

j\k) = ^^p-e,{k)e''"/^°, 

ITiq k" 

where pis the difference between the spatial momenta of a quark and antiquark with 
the mass rUq, t\{k) is the vector of linear polarization of a photon; Eq means the 
phenomenological parameter which depends on the initial energy. For two sources 
the transition current is 

2 

n=l 

Index n labels the independent quark-antiquark interactions taking place at dif- 
ferent space-time points |/„ which are considered to be randomly distributed in 
the space-time volume of the (photon) source. We suppose the gamma-quanta 
are created by some random currents (sources) J^(a;) through the Lagrangian den- 
sity Lint = —gJiJ.{x)A^{x). The neutral weak current of (final) charged leptons is 
Jjj^{x) = —gJ{x)'y^l{x) with / : e,fj,,T, for which the problem can be solved exactly. 
Both the currents J^(a;) and JUx) exist in a restricted space-time region and they 
are chaotically and randomly disturbed by external fields (forces). The restricted 
domain is characterized by the internal stochastic scale Lgt, the meaning of which 
is explained in [2,3]. 



By studying BEC of identical particles, it is possible to determine the time scale 
and spatial region over which particles do not have the interactions. Such a surface is 
called as decoupling one. In fact, for an evolving system such as pp collisions, it is not 
really a surface, since at each time there is a spread out surface due to fluctuations in 
the last interactions, and the shape of this surface evolve even in time. The particle 
source is not approximately constant because of energy-momentum conservation 
constraint. 

Actually, the second order distribution function 

normalized to the product A^i(A;i) ■ N2{k2) with one-particle distribution functions 
Ni{ki) = {\J^'{ki)\'^) {i = 1,2) formally defines the probability to find two photons 
with momenta ki and k2 issued at yi and y2- The crossing momenta has to be taken 
into account. 

In this work, we make an attempt to demonstrate that the problem of properties 
of the genuine interactions can be explored using experimental data which can be 
collected at the LHC These data can be analyzed through the compared measures 
of some inclusive distributions and final state correlations. 

One of the aims of this paper is to carry out the proposal for the experimental 
measurements of virtual 7*7* pair correlations. 

This exploration is theoretically supported by the quantum field theory model 
approaches [4-9,2,3] at finite temperature, (QFTfj), where one of the main parame- 
ters is the temperature of the particle (emitting) source under the random external 
forces (fields) infiuence. 

We propose that photons do not strongly interact with (produced) medium: 
they carry information about early stage of reaction. Our pragmatic definition is: 
photons are produced not from hadronic decays. Any source of real gamma-quanta 
produces virtual photons with very low mass. If the momentum of the virtual 
gamma-quantum is sufficiently small, the source strength should be small as well. 
The real gamma-quantum can be measured from the virtual yield which is observed 
as low mass of the spectrum for lepton-antilepton pair. 

The main channels are the two-photon production pp — )■ Higgs — )■ 7*7* — )■ 
2/i~2/i"'", 2e~2e"'", e~e~^n~fi~^, ... in pp collisions. An efficient selection of leptons 
needs to be performed according to the following criteria. First, all leptons were 
required to lie in the pseudorapidity range covered by, e.g., the CMS muon system 
that is, \ri\ < 2.4. Second, the leptons were required to be unlikely charged in pairs. 

The di-lepton channel is especially promising from the experimental point of 
view, since it is expected that the experimental facilities related for LHC will make 



it possible to record muons of energy in the TeV range with a resolution of about a 
few percent and an efficiency close to 100 %. Moreover, this channel is characterized 
by a maximum signal-to-background ratio in the energy region being considered. 

2. A pair of identical bosons with momenta pi and p2 and the mass m produced 
incoherently from an extended source will have an enhanced probability 6*2(^1,^2) = 
-^i2(Pi,P2)/[-^i(Pi) ■ ^^2(^2)] to be measured in terms of differential cross section a, 
where 

a dill ""2 
to be found close in 4-momentum space 5R4 when detected simultaneously, as com- 
pared to if they are detected separately with 

In an experiment, one can account the inclusive density ^2(^1,^2) which describes 
the distribution of two particles in Q (the sub- volume of the phase space) irrespective 
of the presence of any other particles 

, ,11 (Pn2 

P2{P1,P2) 



2! Uevents dpi dp2 ' 

where n2 is the number of particles counted in a phase space domain {pi + dpi,p2 + 
dp2). The multiplicity A^ normalizations stand as 

p{p)dp={N), 
a 

P2{pi,P2) dpi dp2 = {N{N - 612)), 
in 

where (A^) is the averaged number of produced particles. Here, S12 = for different 
particles, while 612 = 1 in case of identical ones (coming from the same event). 

On the other hand, the following relation can be used to retrieve the BEG func- 
tion C2(Q): 

c (o) - ^(^^ ri) 

^2(Q) - j;^^^^Qy (1) 

where N{Q), in general case, is the number of particle pairs (off-shell photons) in 
BEG pattern with 



Q = \ -{Pi - P2)^ ■ {Pi - P2Y = VM2-4m2. (2) 



In definitions ([T]) and ([2]), N'^'^^ is the number of particle pairs witliout BEC and 



p^- = (ujijPi) are four-momenta of produced photons {i = 1, 2); M = \/{pi +^2)^ 
is the invariant mass of the pair of photons. 

An essential problem in two-particle correlations is the estimation of the ref- 
erence distribution N^'^^ {Q) in Eq. ([1]). If there are other correlations beside the 
Bose-Einstein effect, the distribution N'^'^^{Q) should be replaced by a reference 
distribution corresponding to the two-particle distribution in a geometry without 
BEC. Hence, the expression ([T]) represents the ratio between the number of 7*7* 
pairs N{Q) in the real world and the reference sample N^^^ [Q) in the imaginary 
world. Note, that the reference sample can not be directly observed in an experi- 
ment. Different methods are usually applied for the construction of reference samples 
[1] , however all of them have strong restrictions. One of the preferable methods is to 
construct the reference samples directly from data. The 7*7* BEC can be estimated 
for each bin of the photon average transverse momentum p^ = Ipj'^ + PT2I/2 as the 
ratio of the distribution of photon pair invariant relative momenta where both pho- 
tons with transverse momenta pt^ and /% were taken from the same event to the 
same distribution but with the photons of the pairs taken from different events. For 
our aim, for the reference sample N'^'^^{Q), it is suitable to use the pairs 7*7* from 
different (mixed) events. 

It is commonly assumed that the maximum of two-particle BEC function C2{Q) 
is 2 for pi = P2 if no any distortion and final state interactions are taking into 
account. 

In general, the shape of C2{Q) is model dependent. The most simple form of 
Goldhaber-like parameterization for C2{Q) [10,11] is often used for experimental 
data fitting. 

C2iQ) = Co ■ (1 + Ae-«'^') ■ (1 + eQ), (3) 

where Co is the normalization factor, A is the chaoticity strength factor, meaning 
A = 1 for fully incoherent and A = for fully coherent sources; the symbol R is 
often called as the " correlation radius" , and assumed to be spherical in this param- 
eterization. The linear term in ([2]) is supposed to be account within the long-range 
correlations outside the region of BEC. Note that the distribution of bosons can be 
either far from isotropic, usually concentrated in some directions, or almost isotropic, 
and what is important that in both cases the particles are under the random chaotic 
interactions caused by other fields in the thermal bath. In the parameterization ([3]) 
all of these issues are embedded in the random chaoticity parameter A. To advocate 
the formula (^ it is assumed: 
a. incoherent average over particle source where A serves to account for: 



- partial coherence, 

- 77 purity; 

b. spherical Gaussian density of particle emission cell (with radius R); 

c. static source which means no time (energy) dependence. 

However, to enlarge the quantum pattern of particle production process and to avoid 
the static and undistorted character of particle emitter source, we have already sug- 
gested to use the C2{Q) function within QFTj^ accompanying by quantum evolution 
approach in the form [2,3]: 

C2ipi,P2) ^ aN) {1 + Ai(/3) e-'^''- [1 + A2(/3) e^^'^^'] } , (4) 

where exp(— Aps^) is the smearing smooth dimensionless generalized function with 
^p5R = (Pi— P2)^ ^ij.v {pi~'P2Y ■ ^fiu is the nonlocal structure tensor of the space-time 
size and it defines the domain of emitted photons. ^(A^) depends on the multiplicity 
A^ as ^(A^) = {N{N - l))/{Nf. The functions Ai(/3) and A2(/3) are the measures of 
the strength of BEG between two photons: Ai(/3) = 7(0;, /3)/(l + a)^, and the correc- 
tion to the coherence function in the brackets of Eq. (j4]) is A2(/3) = 2a/ \/^{u],j3). 
The function 7(0;, /3) calls for the quantum thermal features of BEG pattern and is 
defined as 

7(c.,/3)=7H = ^(^)^(^.) , n(a;)=n(c.,/3) = ^^^^-^, c. = ^-, (5) 

where n(w, /9) is the mean value of quantum numbers for Bose- Einstein statistics 
particles with the energy uj in the thermal bath with statistical equilibrium at the 
temperature T = 1/(3. The following condition '^fUfloj, (3) = N is evident, where 
the discrete index / stands for the one-particle state /. 

The important parameter a(/3) in (^, the measure of chaoticity, summarizes our 
knowledge of other than space-time characteristics of the particle emitting source, 
and it varies from to cxd (see [12] for details). 

In terms of time- like Rq, longitudinal Rl and transverse Rt components of the 
space-time size i?^, the distribution Apjf looks like 

A,^ ^ A,n = {Apy Rl + {Ap^ Rl + {Ap^^f R^. (6) 

i?o in ([6]) is treated as the measure of the particle emission time, or even it represents 
the interaction strength of outgoing particles. 

Hence, we have introduce a new parameter R^, a four-vector, which defines the 
region of nonvanishing particle density with the space-time extension of the particle 
emission source. Formula (jl]) must be understood in the sense that exp(— Apsjj) is a 



distribution that in the hmit R -^ oo strictly becomes a. 6 - function. For practical 
using with ignoring the energy-momentum dependence of a, one has: 



C2(g)^e(iV) l + Ai(/3)e 



-Q2/?2 



l + A2(/3)e 



+Q2i?2/2 



(7) 



The parameter R in formula ([7]) is the measure of the space-time overlap between two 
photons, and the physical meaning of R depends on the fitting of C2(Q)-function. 
R can be defined through the evaluation of the root-mean-squared momentum Qrms 
as: 



QlmsiP) = iQ') 



f,^d\Q\Q' C2{Q,P) 



Io^d\Q\ C2iQ,^)-l 



C2iQ,(3) 



C2{Q.P) 



where R and Qrms{P) are related to each other by means of 

1/2 



R = R{(3) 



1 + 



1 



1 



4a(/3) V 2 



^ rms\P j 



We find the following restricted window a/S/^ < {R ■ QrmsiP)) < VS, where the 
lower bound satisfies to the case a — )■ (no any distortion in the particle production 
domain), while the upper limit is given by the very strong influence of chaotic 
external fields (forces), a — )■ oo. The result is rather stable in the wide range of 
variation of a. 

3. It has been emphasized [2,3] that there are two different scale parameters 
in the model considered here. One of them is the so-called "correlation radius" R 
introduced in Q and also presented in (|1]). In fact, this i?-parameter gives the 
pure size of the particle emission source without the influence of the distortion and 
interaction forces coming from other fields. The other (scale) parameter is the scale 
Lgt of the production particle domain where the stochastic, chaotic distortion due 
to environment (the influence of other fields, forces) is enforced. This stochastic 
scale carries the dependence of the particle mass, the a-coherence degree and what 
is very important - the temperature T-dependence. 

One question arises: how can BEG be used to determine the effective scale Lst 
and, perhaps, the phase transition? We suppose the changes of 7*7* production 
region and effects yielding the dynamical variables and parameters of BEG due to 
in-medium distortion. 



Consider the finite system (expanding or squeezing) with the flow of two gauge 
bosons pairs, eg., B^B^^- pairs. The Hamiltonian is 



Ho = - I d^x 



\VB\^ + m^B^, 



which is asymptotically free in the rest frame of undistorted matter with the fleld 
B^ = (0, B) having the mass m in general case. This Hamiltonian and commutation 
relations can possess the exact symmetry. However, the observed states in real 
physical environment can not be realized in the framework of this symmetry. Hq 
has to be added by the Hamiltonian 



-^flU 



Hdist = ^ I d^xd^yB^{x)5F^{x - y) B,{y), 

which is provided by the distortion due to in- medium effect, in particular, because 
of temperature of the environment. The fleld B^{x) propagates in medium with p 



(momentum) - dependent effective (squeezing) frequency Ujb = Jm? + fP — 6F^{p) 
and, consequently, the mass m^ is related to the undistorted (asymptotic) mass m 



by 17113 = \ rrn? — 6F^{p). Here, 6F^{x) is the non-local formfactor leading to the 

modiflcation of the particle frequency (mass) of undistorted matter. Because of the 
quadratic form of the Hamiltonian Hdist through the asymptotic Bose - operators of 
annihilation (creation) (i{p){a~^{p)) in 



d^P ^^ ^x, 

(2 7r)V^^+P^ A 



B,{x) = / ,^ ,^ ^^^ E e;(p)a(p)e-^^ 



one can use the Bogolyubov transformation 

hf = Ufttf — Vf a^r, bt = -Uj at — v*f a-f 

for new operators of annihilation (creation) bf{bt) in distorted medium. In some 
sense, the last operators correspond to thermalized-distorted quasiparticles. The 
functions Uf and Vf obey the condition \uf\ — \vf\ =1 and can be given in the 
form Uf = coshrif, Vf = sinhr^j, where r]f has to be even function of index /. 
For simplicity, we identify / with momentum p, and for squeezing frequency one 
gets rjp = In u"^ (p) / u'i{p) . Indeed, the Bogolyubov transformation is equivalent to 
squeezing procedure. 

4. In the Higgs-boson rest frame, there is a kinematical conflguration for two 
pairs of flnal lepton momenta {pup'j) and {quq'j): Pi — Pj and qi ~ g[. Thus, the 
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final state mimics a two-body final state, and if the leptons being the electrons or 
even muons, the virtual photons with momenta p and p' are nearly on mass-shell, 
p'^ r-^ p' ~ 4m^, where mi is the lepton mass. For the configuration given above, the 
stochastic scale Lgt has different behavior depending on T. At lower temperatures 
we have 



^st 



^3/2 ^2^ ^2m,/T 



A^a{N)mf'^T^/^ 



1 + ^^, 

16 mi J _ 



where the condition 2 ra /3 m/ > 1 is taken into account for any integer n; Mh is the 
mass of the Higgs-boson, m,i stands for the lepton mass. 

On the other hand, at higher temperatures, when T > 2nini, one has 



■'st 



TT^Ml 



_48C{3)a{N)T^mf_ 



C(3) = 5;]n-3 = 1.202. 



(9) 



n=l 



In case of two real photons correlation one can use the longitudinal stochastic 
scale instead of (jHD and (jH]), L°^^{mi = mTJ'i), with the average transverse mass 

my = 0.5(wp^^ -|- \ Pt^) in the frame of, e.g., the Longitudinal Center of Mass 

System (LCMS) [13]. 

It turns out that the scale L^t defines the range of stochastic forces. This effect is 
given by Q;(A^)-coherence degree which can be estimated from the experiment within 
the function C2{Q) as Q close to zero, 6*2(0), at fixed value of mean multiplicity (A^): 



c^iN) = ^—-^ '- , (10) 

where 7(n) is defined in (E]) and (52(0) = C2(0)/^(A^). The upper limit of 6*2(0) 
depends on (A^) and the quantum thermal features of BEC pattern given by 7(n): 
6*2(0) < ^(A^)[l — 7^/^(n)/2]~^. This upper limit is restricted by the maximal value 
of 2 in the ideal case as (A^) — )■ 00 and 7(n) = 1. 

Note, that for C*2((5) - function (jl]), the limit a — ;■ 00 yields for fully coherent 
sources with small (A^), while a — )■ case stands for fully chaotic (incoherent) 
sources as (A^) — ?■ 00. Actually, the increasing of T leads to squeezing of the domain 
of stochastic forces influence, and LgtiT = Tq) = i? at some effective temperature 
Tq. The higher temperatures, T > Tq, satisfy to more squeezing effect and at the 
critical temperature T^ the scale Lst{T = T^) takes its minimal value. Obviously, 
Tc ~ 0(200 GeV) defines the phase transition where the chiral symmetry restoration 



will occur. Since in this phase all the masses tend to zero and a — > at T > T^, one 
should expect the sharp expansion of the region with Lst{T > T^) — t- oo. 

Using the relation between Rl and rriT obtained from the Heisenberg uncertainty 
relations [14] Rl^itit) = Ca//i-Ro/2 vrrrir one can estimate Rq within the formula: 

(11) 



Rq 



c^ h 



Tl^Ml 



_3C(3)a(iV)To3m^ 



where c and h are the speed of light in vacuum and the Planck constant, respectively. 
The expression fITT]) relates the two-photon emission time Rq with the Higgs-boson 
mass Mhi the transverse mass of two photons uit and the temperature of the BEC 
domain Tq > nrriT- From fITT]) one obtains Rq ~ 10~^^sec and Rq ^ 0.45 ■ 10~^^sec 
for To = 50 GeV and Tq = 100 GeV, respectively, at Mh = 120 GeV and tut = I 
GeV in the case of two final muon pairs. The confirmation of this should come from 
the measurement of the Higgs-boson lifetime in two-photon decay channel, as well 
as from the BEC data of electroweak interacting particles. 

The qualitative relation between R and L^t above mentioned is the only one 
we can emphasize in order to explain the mass and temperature dependencies of 
the source size. The dependence of the stochastic measure of chaoticity a on the 
minimum scale cut Lgt and Qrms can be used to define the fit region for different 
Pt- Such a minimum cut on Lgt introduces a lower cutoff on R. 

There are a number of effects which may give rise to Q — correlations and, thus, 
mimic the two-photons BEC coming from Higgs decays. These include i) correlations 
of light hadrons and/or vector bosons misidentified as photons, ii) radiative decays 
of resonances in both pseudoscalar and vector sectors, in) collective flow, etc. 

Apparatus or analysis effects which may result in the sensitivity of external 
random forces influence, may be investigated by studying the dependence of the 
correlation functions on the stochastic scale Lgt- This effect is expected to contribute 
strongly at small Lgt (or large a and T). 

Our model is consistent with the idea of the unification of weak and electromag- 
netic interactions at T > Tc, predicted by Kirzhnitz and Linde in 1972 [15,16]. In 
addition, there is the analogy with the asymptotic free theory: the properties of 
environment (media) are the same as those composed of free particles in the infinite 
volume (Universe). 

5. To summarize: one of the main reasons to study BEC at finite temperature 
is the possibility to determine the precision with which the source size parameter 
and the strength chaoticity parameter(s) can be measured at particle colliders. Such 
investigations provide an opportunity for probing the temperature of the particle 
production source and the details of the external forces chaotic influence. Moreover, 
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one can predict the mass of the Higgs-boson. No systematic theoretical treatment 
of the Bose-Einstein effect within the Higgs decays in 7* - pair production has been 
given so far. 

In this paper, we faced to the model C2{Q, (3) - function in which the contribution 
of N, T, a are presented. This differs from the methods used in, e.g., LEP experi- 
ments based on the approach ([3]). In fact, the latter resembles the traditional way of 
BEC study, however any qualitative interpretations of R, A, e have no been clarified. 
The model proposed in this paper is expected to be sensitive to the temperature of 
the environment and to the external distortion effects. 

We find that the stochastic scale Lgt decreases with increasing T slowly at low 
temperatures, and it decreases rather abruptly when the critical temperature is 
approached. 

We obtain the dependence of the correlation strength functions Ai and A2 on the 
distance Lgt and the maximal value of C2 at minimum invariant momentum Q and 
T. 

We predict for the first time the spatial size of the source (the correlation radius 
R) of two photons originated from the Higgs-boson decay in restricted domain at 
the proper temperature Tq. The dependence of the Higgs-boson mass, the mean 
multiplicity (A^) and the lepton mass is obtained in the form 

^J/5g2m,/5To 

R ' 



for low values of Tq < nm^, while for higher temperatures one has 



R- ^-' 



The correlation radius R increases with heavier Higgs-bosons at large gamma- 
quantum mean multiplicity (A^). Actually, the experimental measuring of R (in 
fm) can provide the precise estimation of the effective temperature Tq which is one 
of the main thermal characters in the 7*7* pair emitter source in the proper leptonic 



decaying channel 7*7* — > //// with the final lepton energy y kf + mf at given a(A^) 
fixed by C2{Q = 0) and (A^). To is the true temperature in the region of multipar- 
ticle production with dimension R = Lst, because at this temperature it is exactly 
the creation of two 7 quanta occurred in decay of Higgs-boson, and these particles 
obey the criterion of BEC. 

The parameter a can be extracted from the experimental data on the two-photon 
7*7* BEC using C2{Q) with Q being close to zero (see ( ITU]) ). This allows one 
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to estimate the Higgs-boson mass, M^, at the temperature Tq. As a quahtative 
illustration we present here an estimation of Mh assuming To = 50, 100 GeV, a = 10 
% and nil = f^ti = 105 MeV. The result is rather sensitive to the correlation radius 
R and Tq: Mh = 167 GeV at Tq = 50 GeV, Mh = 477 GeV at Tq = 100 GeV for 
R = 1.0 fm; Mh = 68 GeV at Tq = 50 GeV, Mh = 194 GeV at Tq = 100 GeV for 
R = 0.7 fm; 

Because of the fact that {^^fi^) pairs originated from 7* decay overlap in space 
and are created in time almost simultaneously, it is natural to expect that there are 
correlation between (/i"*"/!") pairs coming from different 7*'s due to Bose-Einstein 
interference. These effects may also affect the accuracy with which the (/i"'"/i~) pair 
mass can be measured at the LHC. 
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